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The two-dimensional random gauge XY model, where the quenched random variables are mag- 
netic bond angles uniformly distributed within [— rn, rn] (0 < r < 1), is studied via Monte Carlo 
simulations. We investigate the phase diagram in the plane of the temperature T and the disorder 
strength r, and infer, in contrast to a prevailing conclusion in many earlier studies, that the system 
is superconducting at any disorder strength r for sufficiently low T. It is also argued that the super- 
conducting to normal transition has different nature at weak disorder and strong disorder: termed 
Kosterlitz-Thouless (KT) type and non-KT type, respectively. The results are compared to earlier 
works. 
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I. INTRODUCTION 

The XY gauge glass mode! 1 has attracted much inter- 
est in connection to the vortex glass phase of high-T c su- 
perconductors^. In three dimensions (3D) there is a gen- 
eral consensus that the XY gauge glass model exhibits 
a finite-temperature glass transitioni* 3 *^. However, in 
2D there exist conflicting evidences: On the one hand, 
equilibrium studies of defect energy2iiii£i£, Monte Carlo 
(MC) simulations of the root-mean-square current, and 
resistance calculation^ have suggested that no finite-T 
ordering exists. The glass order parameter has further- 
more been analytically shown to vanish at any finite T— . 
On the other hand, the MC studies of the glass suscep- 
tibilitjiii, as well as dynamical simulations of the resis- 
tanc o 12 ' 13 and non-equilibrium relaxation™, have indi- 
cated a possibility of a finite-T transition. Also the MC 
simulations of the helicity modulus in Ref.™ were inter- 
preted as being compatible with a finite-T transition. 

In this paper we study a generalization of the 2D XY 
gauge glass model — the random gauge XY model — where 
both the temperature T and disorder strength r can be 
varied. When r has the maximum value 1, it corresponds 
to the usual XY gauge glass model, while in the opposite 
limit of r = 0, the standard XY model without disor- 
der is recoveredi^iiSiil. It has been proposed that there 
is a Kosterlitz-Thouless (KT) like transition at a finite 
temperature T c when the disorder strength is sufficiently 
small, and that as r is increased T c becomes smaller un- 
til it vanishes as r reaches the critical disorder strength 
r c (< 1) and T c = for r > r c 16 - 17 . However, even if 
the glass order parameter is zero and even if there is no 
finite-T KT transition for r > r c , the existence of a finite- 
T transition with a different character cannot a priori be 
ruled out. 

In the present paper we perform extensive MC simula- 




FIG. 1: Sketch of the phase diagram in the disorder strength 
r and the temperature T (in units of J/fcs with the coupling 
strength J) plane. There exist two superconducting phases at 
low temperatures: the low-T KT phase marked by "SI" and 
the distinct phase "SII" for small r and large r, respectively. 
The high temperature normal phase is marked as "N". The 
solid phase boundary represents the boundary of the low-T 
KT phase whereas the dashed represents the normal to su- 
perconducting transition for larger r. 



tions to study the phase transition of the random gauge 
XY model. It is found that the system is superconduct- 
ing at any r and that the transition from normal to super- 
conducting phase is consistent with a Kosterlitz-Thouless 
(KT) type at weak disorder and a non-KT type at strong 
disorder. In addition to this we suggest that there ex- 
ist two different superconducting phases at sufficiently 
low temperatures separated by a non-KT phase transi- 
tion (see Fig. QJ. For the special case of r = 1, which 
corresponds to the 2D XY gauge glass model, we also 
compute the root-mean-square current in the same way 
as in Ref. 8 and, in contrast to Ref£, again consistently 
find a finite-temperature transition. 
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II. MODEL AND SIMULATIONS 

The Hamiltonian of the 2D random gauge XY model 
on an L x L square lattice under the fluctuating twist 
boundary condition (FTBC)i£ is given by 



H = — cos 



f 

-I" 



A — Ai 



(1) 



where J is the coupling strength (set to unity from now 
on), the sum is over nearest neighbor pairs, ry = x 
(y) if 3 = i + x(y)- The phase angle Qi at the lattice 
point i satisfies the periodicity Oi+Lx = @i+Ly = @i, an d 
Aij e r[— 7r, 7r] is a uniform quenched random variable 
with the disorder strength < r < 1. The twist variable 
A = (A x , A y ) corresponds to the global twist across the 
system, i.e., the summation of the gauge invariant phase 
difference tfiij along the x (y) direction equals A x (A y ). 
For a given disorder realization, we first compute the dis- 
tribution -P(A), which is related to the free energy F by 
dF/dA = -T(91nP/0A). The twist variable A which 
minimizes F (or maximizes P) is determined from P and 
then fixed when the helicity modulus T = d 2 F/dA 2 and 
the 4th order modulus T4 = d 4 F/dA 4 are computed. 

To ensure that the cooling is slow enough, we simul- 
taneously cool two replicas (a and (3) of the system and 
measure Aq'' 3 . For the first cooling at a new temper- 
ature we use f 20000 update sweeps (for spin and twist 
variables respectively). Then we check that 
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where 5 sets the precision of the cooling. The idea with 
the annealing condition Eq. J2J is to keep the system close 
to the lowest-energy state at a particular temperature (so 
the system does not freeze into a local minimum that bias 
[T] and [T4]). Since the system moves more swiftly over 
the configuration space the higher the temperature is we 
can choose 6 increasing with temperature. A choice that 
proves good in practice is 




for T < 0.3 

for 0.3 < T < 0.6 

for T > 0.6 



(3) 



If the annealing condition Eq. @ fail we repeat the cool- 
ing with three times as many update sweeps; if it fails 
again we increase the number of cooling sweeps a fac- 
tor three again, and so on until the condition is fulfilled. 
When Ao is chosen, before cooling, we let the system run 
until Eq. J5J) is fulfilled with A replaced by A. 

We repeat the above calculations for more than 500 
different disorder realizations (the disorder average is de- 
noted by [• • •] throughout this paper) . 

For the case of r — 1, we also use the periodic boundary 
condition PBC, corresponding to A = in Eq. Q, and 
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FIG. 2: (a) The size dependence of the helicity modulus 
[T] at different temperatures T for the fully disordered case 
(r — 1) — the gauge glass model. At high T, [T] is shown to 
vanish as L is increased, whereas it saturates to a finite value 
at low T (qualitatively given by the dashed lines are linear ex- 
trapolations from the smallest sizes and correspond to power 
law behavior) (b) Finite-size scaling of [T] of the data in (a). 
From the standard finite-size scaling form, the critical tem- 
perature Tc ~ 0.19 and the critical exponent v ~ 1.1 are 
determined. 



compute the root-mean-square current defined by& 



Jrms — 




sm( 



1/2 



(4) 



for comparisons with earlier works. 

For J rms , for each disorder, we use 50000 update sweeps 
for thermalizations and 500000 sweeps for measurements 
where the actual measurement was performed every tenth 
sweep. We used this computer time saving strategy since 
we found that performing the measurement every sweep 
gave closely the same result. 
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FIG. 3: The root-mean-square current / rms . (a) The finite- 
size-scaling plot used in Ref, 8 for T c = shows systematic 
deviations in the low-temperature region (y = 2.2 in Ref. 8 
is used). We include the data in Ref£ for comparisons, (b) 
Finite-size scaling form in Eq. @ yields T c » 0.2 and b « 0.5. 
Inset: All data points in the main panel collapse to a smooth 
curve with v = 1.1, the same value as found for the helicity 
modulus in Fig. [5] 



III. RESULTS FROM SIMULATIONS 



We first investigate the standard 2D XY gauge glass 
model, corresponding to the fully disordered case (r = 1), 
and show in Fig. EJa) the helicity modulus [T] as a func- 
tion of the system size L for various temperatures. It 
is clearly shown that at high T, [T] goes to zero as the 
system size L is increased. The crucial point is that, at 
low enough temperatures (T < 0.2), [T] changes its cur- 
vature in terms of L, and appears to saturate to a finite 
value as L is increased. This behavior suggests a phase 
transition with a scale-invariant power law dependence 
at the critical temperature and a diverging characteristic 
length. Such a behavior can often be described by the 
finite-size scaling form 



where [T] ~ L~ b at the critical temperature T c , and the 
critical exponent v is related to the divergence of the 
coherence length. At T c , the scaling function / has the 
same value irrespective of L, implying that L b [T] versus 
T should have a unique crossing point at T c for various 
sizes. In Fig.^b), it is clearly shown that [T] has a scale- 
invariant behavior [T] oc L~ at a unique T, signaling 
a phase transition. The inset of Fig. |2fb) furthermore 
confirms that this scaling behavior is consistent with the 
standard form JSJ with v « 1.1. 

Simple dimensional analysis for the non-disordered 
case, r = 0, gives for [T] the exponent 6 = and from 
such a dimensional analysis one would likewise conclude 
that 6 = also for the disordered case. Thus if [T] scales 
with 6 ^ for the disordered case this is equivalent to the 
appearance of an anomalous dimension not accounted for 
by simple dimensional analysis. Our suggestion, based on 
the simulation results, is consequently that the disorder 
introduces such an anomalous dimension. 

Figure [3] shows the root-mean-square current Q for 
the PBC. We first note in Fig. EJa) that the finite- 
size-scaling form for T c = used in Ref£, J rms = 
L-i/v/fTL 1 /") with v = 2.2, shows systematic devia- 
tions from the data collapse to a single scaling curve at 
lower temperatures, in contrast to what was concluded in 
Ref A, when more and better converged data are included. 
Furthermore a T c — 0-collapse cannot be achieved with 
any value of v. On the other hand if we, in analogy with 
the finding for [T] above, use the scaling form 



>ff L V»(T-T c ) 



(0) 



[T] = L~ h f (j J 1 / u (T — T c 



(5) 



which allows for the same anomalous dimension, one ob- 
tains the scaling plot in Fig. OJb) with 6 « 0.5 and 
T c « 0.2i& From Figs. [2 and |3 we conclude that the 2D 
XY gauge glass (r = 1) exhibits a non-KT type super- 
conducting to normal transition at T c w 0.2 characterized 
by the existence of the anomalous dimension 6 sa 0.5 and 
v w 1.1. 

Next we investigate the phase transition with decreas- 
ing r. For r — 0.9, 0.8, • • • , 0.5, we obtain the finite-size 
scaling plots of the same quality as in Fig. OJb) , and de- 
termine the phase boundary (the dashed line in Fig. 
As r is changed from 0.5 to 0.4 the exponents 6 and v 
exhibit quite rapid changes, 6 from 0.27 to 0.06 and v 
from 1.1 to 2.0. This, in our interpretation, reflects that 
near r = 0.4 the nature of the superconducting-normal 
transition changes from a non-KT type to a KT type. 

For r = the phase transition is of the KT nature and 
it has been suggested that this character should persist 
along the phase boundary up to some r c < l 16,17 . The 
KT transition is characterized by that [T] jumps from 
a finite value [T] = 2T/ir to zero as the phase line is 
crossed from the small (T,r) region. For r — we find 
that the KT transition is also characterized by the in- 
crease of the 4th order modulus |[T 4 ]| oc L c at T c with 
a positive exponent c, and [T4] stays at a constant value 
below T c as L is increased^. Figure 01 (a) shows [T4] 
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FIG. 4: The fourth order modulus and helicity modulus (in- 
set) as a functions of r for temperatures (a) 0.6 and (b) 0.1. 
The solid line in the inset represents the universal-jump con- 
dition for a KT transition [T] = 2T/tt. 



for T = 0.6 as a function of r. The KT condition in the 
inset gives r ~ 0.4. This means that if the transition is 
of KT type then it would occur around r w 0.4. The on- 
set of size dependence for [T4] in Fig. 0] (a) is consistent 
with a transition at around r ^ 0.4. Thus the transition 
is compatible with a KT-character, although a different 
character cannot be ruled out. The situation for T = 0.1 
in Fig. H|(b) is quite different: the inset shows that the 
KT jump condition is not fulfilled and a KT transition 
can be ruled out. Yet there is a marked structure in [T] 
around r ~ 0.4. This structure corresponds to the onset 
of strong size dependence in [T4] . We interpret this onset 
as the reflection of a true divergence in [T4] consistent 
with a phase transition. Thus we suggest that the whole 
boundary to the low-T KT phase (solid line in Fig. ^| is 
reflected in a divergence of [T4] and that this line ends 
at (T = 0,r c w 0.4). A phase line which ends at (T = 0, 
r c « 0.4) has been found in many earlier investigations 
(see e.g. Refii& and references therein)Si. The differ- 
ence with earlier work is that in our case such a phase 
line is for lower temperatures between the two different 
superconducting phases SI and SII (as demonstrated by 
our direct calculation of [T]), whereas earlier work have 



concluded that such a phases line does indeed exist but 
separates a superconducting phase from a normal phase 
all the way down to T = 0. Our suggested phase diagram 
is thus consistent with earlier work as to the existence of 
a phase line ending at (T — 0, r c 0.4). For tempera- 
tures larger than the merging with the second phase line 
(above dashed line in Fig. the transition is consistent 
with a KT transition although a different character can- 
not be ruled out. However, below the merging with this 
second line the transition is not a KT transition. It may 
be that there still is a jump in [T] at this transition, but 
this is then between two non- vanishing values. 

Based on the numerical evidences we suggest the struc- 
ture of the phase diagram is sketched in Fig.^ One strik- 
ing feature is the finite-T transition line between normal 
and superconducting phases, starting from the boundary 
of the low-T KT-phase (solid line) and changing into the 
(dashed) line which ends at T c « 0.2 at r — 1. The lat- 
ter line (dashed line) is characterized by the appearance 
of an anomalous dimension b. Since the KT transition 
does not have such an anomalous dimension it follows 
that, if the transition along the boundary to the low-T 
KT phase boundary has KT character, then b should ap- 
proach 6 = when the two transition line merge. Thus 
the rapid drop from b w 0.27 to b w 0.06 as r is changed 
from r = 0.5 to 0.4 is consistent with a change over to 
a KT transition. Another interesting feature is that the 
phase line (solid line in Fig-HJ, associated with the diver- 
gence of [T4], continues even inside the superconducting 
region and ends at r « 0.4 for T = 0. Although the tran- 
sition separating the KT phase and the normal phase (SI 
and N, respectively, in Fig. [JJ, may well be of the true 
KT type, manifested by a universal jump in the helicity 
modulus, the phase line separating SI and SII in Fig. 
is not a KT transition and the helicity modulus has a 
non-zero value on both sides. 

It is interesting to note that earlier works have found 
evidences for only two phases separated by a single phase 
line; either, in the more prevailing view, a phase line end- 
ing at a point T = for a finite r c < 1, or, in the less 
prevailing view, a phase line ending at a point T > for 
r = 1. ^From our numerical simulations we instead sug- 
gest three distinct phases separated by two phase lines, 
which combines the two earlier proposed scenarios and 
provides a unified picture: On the one hand, in Ref. 5 , 
the end point of the phase boundary to the low-T KT- 
phase has been obtained to be r c « 0.37 and T = 0, 
which is consistent with our cruder estimate r » 0.4. On 
the other hand, the phase transition point at T ks 0.22 
with v — 1.1, found in Refs^i ancU^ for r = 1, is in very 
good agreement with the end point of our finite-T phase 
line (T » 0.2, v w 1.1). However, the difference with the 
previous work 5 is that according to our interpretation the 
phase line below T w 0.2 ending at T = and r c rj 0.4 
separates two distinct superconducting phases, while the 
whole phase line in Ref* 5 " is for the superconducting to 
normal transition. A very hand-waving picture of the 
scenario of our phase diagram in terms of vortex motion 
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is sketched in Fig. In SI the vortex motion is sup- 
pressed by vortex pair binding, in N pair unbinding has 
occurred and free vortices exists which are not entirely 
pinned by the disorder, whereas in SII vortex pair un- 
binding has occurred but the vortices are pinned by the 
disorder. 



IV. SUMMARY 

One main conclusion from this phase diagram is that 
the XY gauge glass model (r = 1) has a finite-T tran- 
sition. How is this possible in view of earlier conflicting 
evidences? In Ref£ it was concluded, on the basis of 
an analysis of data for the root-mean-square current for 
standard periodic boundary conditions 2 ^, that no finite- 
T transition exists in 2D. We have found that, taking into 
account the possibility of an anomalous scaling dimen- 
sion, / rms displays a transition at a finite-temperature 
(see Fig.|3J^. Another puzzling evidence to the contrary 



is the T = calculations of the size scaling of the defect 
energyiiSii. However, as discussed in Refsi and^ii, the 
local vorticity conservation must be properly taken into 
account when calculating the energy barriers. Thus the 
energy barrier for vortex dissipation increases with sys- 
tem size when taking the local vorticity conservation into 
account^. This growing of the energy barrier for vortex 
dissipation with system size supports the possibility of a 
finite-T transition 23 . 

The appearance of a new superconducting phase for 
the XY random gauge model is intriguing. In particular, 
since it is neither a low T KT phase nor a phase with 
a finite glass order parameter. The true nature of this 
phase and the existence of similar phases in other related 
models are open questions. 
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